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Abstract. The problem of optimal control of non-Markovian open quantum system 
via weak measurement is presented. Based on the non-Markovian master equation, 
we evaluate exactly the non-Markovian effect on the dynamics of the system of 
interest interacting with a dissipativc reservoir. We find that the non-Markovian 
reservoir has dual effects on the system: dissipation and backaction. The dissipation 
exhausts the coherence of the quantum system, whereas the backaction revives it. 
Moreover, we design the control Hamiltonian with the control laws attained by 
the stochastic optimal control and the corresponding optimal principle. At last, 
we considered the exact decoherence dynamics of a qubit in a dissipativc reservoir 
composed of harmonic oscillators, and demonstrated the effectiveness of our optimal 
control strategy. Simulation results showed that the coherence will completely lost in 
the absence of control neither in non-Markovian nor Markovian system. However, the 
optimal feedback control steers it to a stationary stochastic process which fluctuates 
around the target. In this case the decoherence can be controlled effectively, which 
indicates that the engineered artificial reservoirs with optimal feedback control may 
be designed to protect the quantum coherence in quantum information and quantum 
computation. 
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1. Introduction 

Quantum information science has emerged as one of the most exciting scientific 
developments of the past decade. The most difficult problem in realizing the quantum 
information technology is that the quantum system can never be isolated from the 
surrounding environment completely. Interactions with the environment deteriorate the 
purity of quantum states. This general phenomenon, known as decoherence [H [2] , is a 
serious obstacle against the preservation of quantum superpositions over long periods of 
time. Decoherence entails nonunitary evolutions, with serious consequences, like a loss 
of information and probable leakage toward the environment. Thus, on the one hand, 
the information carried by a quantum system has to be protected from decoherence. 
On the other hand, a detailed knowledge of quantum dynamics and control will help 
to construct high precision devices that accurately perform their intended tasks despite 
large disturbances from the environment. 

For the purpose of efficiently processing quantum information, results and methods 
coming from the classical control theory have been systematically introduced to 
manipulate practical quantum systems [SI HI [H El El El [HI [TOl El US [13] . Early works aim 
to answer the questions like controllability and observability in finite dimensional closed 
quantum systems p31 [15] . They have many applications like molecular dynamics in 
quantum chemistry. As we demonstrated above, quantum systems are very vulnerable 
when exposed in a noisy environment. So, open quantum systems control is becoming a 
matter of concern for more and more physicists and cyberneticists. Recently, quantum 
control works include: (i) steering a quantum system from its initial state to a given 
final state or a set of final states [16]. Transferring to specific final states are very 
important for applications to quantum computing, quantum chemistry and atomic 
physics, (ii) Quantum decoherence control [3l [11]. Suppression of decoherence is 
essential to effectively protect quantum purity and quantum coherence, (iii) Quantum 
entanglement control [171 [ISj. Reliable generation and distribution of entanglement in 
a quantum network is a central subject in quantum information technology, especially 
in quantum communication. 

To the quantum control strategies, a number of interesting schemes have been 
proposed during the last few years in order to protect the purity and counter the effects 
of decoherence. It is just like the classical control theory, according to the principle 
of controllers' design, the open quantum systems control includes open-loop control 
and closed-loop control. Quantum open-loop control means to design the controller 
without measurement. Early proposed quantum error-correction codes, error- avoiding 
codes and Bang-Bang control can be classified as open-loop control. Nowadays, by 
tuning the system's Hamiltonian, coherent control theory has opened new perspective 
on decoherence and entanglement control. It can decouple a part of the system 
dynamics from decoherence, on which noiseless quantum information can be encoded 
either in Markovian open quantum systems or in non-Markovian open quantum system 
[T9l [20] . The quantum coherent control plays an important role in state generation 
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and transfer in quantum information, quantum chemistry and optical physical, and 
it is a focus these days. Feedback is the essential concept of classical control theory, 
but the study of feedback control has only a toehold in quantum control theory for 
many fundamental problems needed to be solved. The main problem is quantum 
measurement. Feedback control needs measurement to reduce the system's uncertainty, 
whereas measuring a quantum system will inevitably lead to quantum state collapse, 
which increase the system's uncertainty in turn [21]. As we demonstrated above, the 
theory of quantum measurement is strange in that it does not allow the noncommuting 
observables to be measured simultaneously. The essence of feedback control is expected 
to attain robustness to noise or modeling error, and quantum control using continuous 
measurement, so-called measurement based quantum feedback control, was proposed 
in early 90's. Recently, the stochastic control theory has been exploited to the open 
quantum system. The theory of quantum feedback control with continuous measurement 
can be developed simply by replacing each ingredient of stochastic control theory by its 
noncommutative counterpart. The system and observations are described by quantum 
stochastic differential equations. The field of quantum stochastic control was pioneered 
by V. P. Belavkin in the remarkable paper [22] in which the quantum counterparts of 
nonlinear filtering and LQG control were developed. The advantage of the quantum 
stochastic approach is that the details of quantum probability and measurement are 
hidden in a quantum filtering equation. In view of physicists, they also reformulated 
the evolution of damped systems with output in the form of an explicitly stochastic 
evolution equation, which specifies the quantum trajectory of the systems. Feedback 
control of quantum mechanical systems which take into account the probabilistic nature 
of quantum measurement is one of the central problem in the control of such systems 
in both physics and control theory. Recently, superconducting qubits [23l |23] have 
also been studied as ways to control and interact with naturally formed quantum two- 
level systems in superconducting circuits. The two-level systems naturally occurring in 
Josephson junctions constitute a major obstacle for the operation of superconducting 
phase qubits. Since these two-level systems can possess remarkably long decoherence 
times. References [25l|26l|27] showed that such two-level systems can themselves be used 
as qubits, allowing for a well controlled initialization, universal sets of quantum gates, 
and readout. Thus, a single current-biased Josephson junction can be considered as a 
multi-qubit register. It can be coupled to other junctions to allow the application of 
quantum gates to an arbitrary pair of qubits in the system. These results [251 EEl ET] 
indicate an alternative way to control qubits coupled to naturally formed quantum two- 
level systems, for improved superconducting quantum information processing. Indeed, 
these predictions have been found experimentally in [28]. More recently, reference [29] 
applies quantum control techniques to control a large spin chain by only acting on 
two qubits at one of its ends, thereby implementing universal quantum computation 
by a combination of quantum gates on the latter and swap operations across the 
chain. They [29] show that the control sequences can be computed and implemented 
efficiently. Moreover, they discuss the application of these ideas to physical systems 
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such as superconducting qubits in which full control of long chains is challenging. 

The rest of the paper is organized as follows. In Section II, we introduce the open 
quantum system and the non-Markovian quantum master equation, and present the 
main difference between the non-Markovian quantum system and the Markovian one. 
The optimal control problem via quantum measurement feedback is precisely formulated 
in Section III. In Section IV, we specifically investigated the stochastic optimal control 
of Spin-Boson system. A useful numerical example is demonstrated in Section V, and a 
few concluding remarks are given in Section VI. 

2. Open quantum system and non-Markovian master equation 

The theory of open quantum system describes the dynamics of a system of interest 
interacting with its surrounding environment [1], and the quantum master equation 
governs the evolution of the quantum system, which plays an important role in relaxation 
and decoherence theory. Markovian approximation is usually used in this master 
equation under the assumption that the correlation time between the systems and 
environments is infinitely short. However, in some cases an exactly analytic description 
of the open quantum system dynamic is needed. Especially in high-speed communication 
the characteristic time scales become comparable with the reservoir correlation time, and 
in solid state devices memory effects are typically non negligible. So it is necessary to 
extensively study the non-Markovian case. In this paper, we will consider the optimal 
control of non-Markovian open quantum system via feedback, typically the quantum 
decoherence optimal control. 

2.1. Open quantum system 

Consider a quantum system S embedded in a dissipative environment B and interacting 
with a time-dependent external field, i.e., the control field. The total Hamiltonian has 
the general form 



(either open-loop or closed- loop), and Hb the bath and Hi the interaction between the 
free system and the bath. Note that the interaction between the control field and the 
free system or the bath are always negligible. The operator Hc{t) contains a time- 
dependent external field to adjust the quantum evolution of the system. One of the 
central goals of the theoretical treatment is then the analysis of the dynamical behavior 
of the populations and coherences, which are given by the elements of the reduced 
density matrix, defined as 



where ptot is the total density matrix for both the system and the environment. For 
simplicity, it is always assumed as the system of a single atom which configured such that 




Ps{t) = tl B[ptot{t)], 



(2) 
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a transition between only two states is possible. The Hamiltonian of the environment is 
assumed to be composed of harmonic oscillators with natural frequencies Ui and masses 

where (xi, X2, - ■ ■ ,xn,Pi,P2, ■ ■ ■ ,Pn) are the coordinates and their conjugate momenta, 
and the Planck constant h = 1 in atomic units and the initial state of the system is 
p(0) = po (for simplicity we write p as ps)- For convenience, we always assume that the 
evolution starts from to = 0. The interaction Hamiltonian between the system S and 
the environment B is assumed to be bilinear, 

Hint = aJ2'^n^Bn. (4) 

n 

The interaction Hamiltonian in the interaction picture therefore takes the form 

= aZnMt)^Bn{t), ^ ^ 

where 

The effect of the environment on the dynamics of the system can be seen as a interplay 
between the dissipation and fluctuation phenomena. And it is the general environment 
that makes the quantum system lose coherence (decoherence). The system-environment 
interaction leads to non-unitary reduced system dynamics. 



2.2. Quantum master equations 

The analysis of the time evolution of open quantum system plays an important role 
in many applications of modern physics. With the Born-Markovian approximation the 
dynamics is governed by a master equation of the relatively simple form 

|p(t) = -t[Hsit),pit)] + ^7™I?[C™]p(t), (6) 

m 

with a time-independent generator in the Lindblad form, her the superoperator 'DlLjp = 
LpL^ — ^L^Lp — ^pL'^L. This is the most general form for the generator of a quantum 
dynamical semigroup. The Hamiltonian Hs{t) describes the coherent part of the time 
evolution. Non-negative quantities 7^ play the role of relaxation rates for the different 
decay modes of the open system. The operators Cm are usually referred to as Lindblad 
operators which represent the various decay modes, and the corresponding density 
matrix equation ([6]) is called the Lindblad master equation. The solution of Eq. ([6]) 
can be written in terms of a linear map V{t) = exp(£t) that transforms the initial state 
p(0) into the state p(t) = V(t)p{0) at time t. The physical interpretation of this map 
V(t) requires that it preserves the trace and the positivity of the density matrix p(t). 
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The most important physical assumption which underhes the Eq. (E]) is the vahdity 
of the Markovian approximation of short environmental correlation times [17] . With this 
approximation, the environment acts as a sink for the system information. Due to the 
system-reservoir interaction, the system of interest loses information on its state into 
the environment, and this lost information does not play any further role in the system 
dynamics. However, if the environment has a non-trivial structure, then the seemingly 
lost information can return to the system at a later time leading to non-Markovian 
dynamics with memory. This memory effect is the essence of non-Markovian dynamics 
[3T| [32| [33] , which is characterized by pronounced memory effects, finite revival times 
and non-exponential relaxation and decoherence. Non-Markovian dynamics plays an 
important role in many fields of physics, such as quantum optics, quantum information, 
quantum chemistry process, especially in sohd state physics. As a consequence the 
theoretical treatment of non-Markovian quantum dynamics is extremely demanding. 
However, in order to take into account quantum memory effect, an integro-differential 
equation is needed which has complex mathematical structure, thus prevent generally 
to solve the dynamics of the system of interest. An appropriate scheme is the time- 
covolutionless (TCL) projection operator technique which leads to a time-local first 
order differential equation for the density matrix. 

The general structure of the TCL master equation is given by 

|p(t) = -z[Hs{t),p{t)] + ^ A„(t)P[C^(t)]p(t). (7) 

m 

The first term describes the unitary part of the evolution. The latter involves a 
summation over the various decay channels labeled by m with corresponding time- 
dependent decay rates Am(t) and arbitrary time-dependent system operators Cmif). 

In the simplest case the rates as well as the Hamiltonian Hs and the operators 
Cm are assumed to be time-independent, that is, it is the Markovian case. Note that, for 
arbitrary time-dependent operators Hsif) and Cm(t), and for Am(t) > the generator 
of the master equation ([7]) is still in Lindblad form for each fixed time t, which may be 
considered as time-dependent quantum Markovian process. However, if one or several of 
the Am(t) become temporarily negative which expresses the presence of strong memory 
effects in the reduced system dynamics, the process is then said to be non-Markovian. 

2. 3. Quantum feedback control 

Quantum feedback control was formally initiated by Belavkin's work [15], [22] and so 
on. The main problem of optimal quantum feedback control is separated into quantum 
filtering which provides optimal estimates of the stochastic quantum operators and then 
an optimal control problem based on the output of the quantum measurement. As we 
demonstrated in the introduction, the quantum system is not directly observable. Since 
the very beginning of the quantum mechanics, the measurement process has been a 
most fundamental issue. The main characteristic feature of the quantum measurement 
is that the measurement changes the dynamical evolution. This is the main difference 
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of the quantum measurement compared to its classical analogue. When one observes 
incompatible quantum events (self-adjoint orthoprojectors with P"^ = P = P* {* denotes 
the Hermitian adjoint) acting in some Hilbert space "H), the state of the system needs to 
be updated to account for the change to the system or back-action. This state change 
was traditionally described by the normalized projection postulate 

which also ensures instantaneous repeatability of the observed event corresponding to 
the projection Pj. 

We now couple the system to a control field. If we assume no scattering between 
the measurement and control fields and assume a weak coupling such that information 
is not lost into the control field, then this effectively replaces the Hamiltonian Hs{t) 
of the system with a controlled Hamiltonian Ho{t) + H{u{t)) for admissible real valued 
control functions u{t) G M. The term feedback refers to control law which is also 
a function of the state p. This control Hamiltonian generates the controlled unitaries 
Ut{u{t)) given the controlled flow Uf{uJ^)pUt{uJJ, where = {u{t)\tQ < t < T} is 
the control process over the interval [0,T). In classical control, we can allow complete 
observability of the controllable system, so that feedback controls are determined by the 
system variables. However, in quantum system the measurement are weak measurement 
(for non-observability of the system operators). So, the feedback controls should be given 
by a function of the stochastic output process lyp^(t), and measurement trajectory is 
fed into the controls Wpiit) '^t(}VF-{'t))- 

Considering both controls on the systems and simultaneous quantum weak 
measurements over multiple observable {Fi}, we have following stochastic master 
equation and measurement outputs for the quantum feedback control system: 

r 

dp = -i[Ho,p]dt + Y,Mt)H[Hk,p])dt + Y,^m{t)V[Cm]pdt 

k=l m 

+ J2 MiV[Fi]pdt + J2 ^/M^n[F,]pdWF^ (8) 

i i 

dWp, = dYi - ^/^,ii{Fip}dt (9) 

where H^-, k = 1,2, ■ ■ ■ ,r are the control Hamiltonians adjusted by the time-dependent 
control laws Uk(t) G M and the superoperator "HfAjp = Ap + pA — [tT{Ap + pA)]p. Note 
that, the last two terms in Eq. ([8]) come from the measurement-induced disturbance 
(we will explain the degenerate in the next section from the quantum measurement 
view of point), where Mj are the effective interaction strengths and r]i are the detection 
efficiency coefficients. Yi are the observation process. In fact, the innovation process 
Wp, is a Wiener process, which satisfy following rules. The expectation of any integral 
over dWp^ or dWp^ vanishes. The differentials dWp^ , dWp. commute with any adapted 
process. Finally, the quantum Ito rules are 

dWp.dW; = dp^^F.dt, dWl. = {dW^f = dW^^.dWp^ = 0. 
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3. Stochastic optimal control of Spin-Boson system 

We have described the quantum feedback controls and have obtained controlled master 
equations, the remaining question is: how do we choose the control strategy u{t) to 
achieve a particular control goal? For a fixed function u{t) the expression ([8]) becomes a 
stochastic differential equation and its solution p*{t) is called a trajectory of the quantum 
system relevant to the quantum control u{t). In general, any trajectory p*{t) should 
fulfill some initial conditions. In this section we approach this problem using optimal 
control theory, where the control goal is expressed as the desire to minimize a certain 
cost functional J{t,u,p). 

The optimality of control is judged by the expected cost associated to the admissible 
control process uj^ for the finite duration T of the experiment. Admissible control 
strategies are defined as those uf^ for which the operator valued cost integral 

J[u{-)]= ! L{T,u,p)dT + 9-g{pT,p*), (10) 

J to 

and the solutions of Eq. ([8]) exist. Here L is a function of time t, the current control law 
u, and the quantum state p. The function g, known as a target or bequest function in 
control theory, is a function of the state parameters at termination. We assume that both 
are continuous in their arguments. 6 is the constant parameters to adjust the penalty 
functional between the two functions. The cost J will vary from one experimental trial 
to another, and must be thought of as a random variable depending on the measurement 
output. 

Here we consider a situation in which a single qubit interacts with a heat bath of 
bosons. We will make use of the two-level non-Markovian master equation. The kinetic 
equation of two-level strong coupling non-Markovian quantum system is the particular 
form of the Eq. ([7]) and ([8]), which performed a rotating wave approximation after 
tracing over the environment without neglecting the counter-rotating terms. By tracing 
out the bath degrees of freedom, we find for p the feedback controlled two-level quantum 
system non-Markovian evolution equation 



2 2 2 

dp= - -Uo[<J^, p]dt - -u^{t)[cT^, p]dt - - 
+ Ti{t)V[a~]pdt + T2{t)V[a+]pdt + MV 

2 J 



Uy{t)[ay, p]dt 
pdt 



' 2 



+ ^/^n [-^1 pdW = fit) 



The free Hamiltonian of the one-qubit quantum system can be written as Hq = 
^(|0)(0| — |1)(1|) = ^CTz, where |0) and |1) are Dirac notations of its two eigenstates 



and Uq is the Rabi frequency, 



ax + la. 



a' 



^[ax - lay] 



with a 



X , ay , (T^ 



the 



Pauh matrices. Here, the controlled part has two channels Hc(t) = ^Ux(t)ax + \uy(t)ay. 
The effective interaction strengths between the measurement and the system is M > 0. 
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Suppose p = \ \ "'""''^ ^^1, thus corresponding to (ITTl) . we have the 
\ X + ly 1 — 2 j 

following evolution equations, 

dx = ^—A(t)x — —X — uoy + Uy{t)z^ dt + \/ MrjxzdW 

dy = (^ojqX - A(t)y - ^y - u^{t)z^ dt + y/WriyzdW 

dz = 2 (^-^x + ^y - Ait)z - 7(t)) dt 

+ ^/M^{-1 + z^)dW. (12) 



Thus the coherence factor A(t) = y^xM-^/2, and the population poo(^) = (1 + 
z)/2, pn{t) = il-z)/2. 

Here, the objective is to compute an appropriate feedback control functions 
Ux(t), Uy{t) steering the system from the initial state po into a target state pr at 
final time T. The relative simple cost functional may be written as 

JHt)] = U\p{T)-pTr + ]- r (uiit) + uirnt, (13) 



2"' ' ' ' " 2 



to 



where || ■ || is the Frobenius norm: = tr = Ylij Here, the first term 

represents the deviation between the state of the system at final time p(T) and the 
target state px, the second integral term penalizes the control field with ^ > 0, the 
weighting factor used to achieve a balance between the tracking precision and the 
control constraints. Minimizing the first term is equivalent to maximizing the state 
transfer fidelity F = tr{p(T)pr}. Our overall task is to find control Uxit), Uy(t) that 
minimizes J[u{t)] and satisfies both dynamic constraint and boundary condition. The 
optimal solution of this problem will be obtained using the optimal principle. The 
corresponding Hamiltonian function may be present in the form 

Hp, u, A) = E |^(«1 + ul) + A^(t)/(t)| , (14) 

where the adjoint state variable X{t) = [Xi{t) , X2(t) , X2,{t)]^ is the Lagrange multipliers 
introduced to implement the constraint. The optimal solution can be solved by the 
following differential equation with two-sided boundary values: 

p = -Q^, p(^o) = Po; (15) 
A=-— , A(T) = p(T) - p^. (16) 



Together with 

m 

du 

mi 



dp 

dn{p*{t),u*{t),x{t),t) 

du 

d'n{p*{t),u*{t),x{t),t) 



= 0; (17) 

<0. (18) 
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Take this into Eq.f ll2p . we get the optimal control for the Spin-Boson system with 
Hamiltonian function Eq.([T] 



U.it) = ^ {\2{t)Xs{t) - X3it)x2{t)} , (19) 



^s/W = o {Ht)Mt) - Ai(t)a:3(t)} . (20) 



e 

2 

The influence of environmental disturbances include both relaxation and dephasing 
effects, which are formulated as Ti(t)V[a~]p and T2(t)V[a^]p. The time dependent 
coefficients Ti^t) = A(t) + 7(t) and r2(t) = A(t) - 7(t), with A{t) and 7(t) being the 
diffusive term and damping term respectively, which can be written up to the second 
order in the system- reservoir coupling constant, as follows, 

A{t) = [ dTk{T)cos{uQT), (21) 
Jo 

7(t) = / dr/i(r)sin(a;or), (22) 

^0 



with 



k{r) = 2 dujJ{uj)coth[nw/2kBT]cos{ujT), (23) 
Jo 

POO 

/x(r) = 2/ dujJ{u)sm{uT), (24) 



^0 

being the noise and the dissipation kernels, respectively. The properties of the Eq. (fTT!) 
strongly depend on the behavior of the dissipation and the noise kernel which, in turn, is 
determined by the spectral density J(cj). In order to obtain true irreversible dynamics 
one introduces a continuous distribution of bath modes and replaces the spectral density 
by a smooth function of the frequency u of the bath modes. In particular we will consider 
the Ohmic spectral density with a Lorentz-Drude cutoff function, 

^(-) = (25) 

where 70 is the frequency- independent damping constant and usually assumed to be 1. 
Uc is the high-frequency cutoff. The analytic expression for the coefficients 7(t)and A(t) 
are given in [TTl [18] . 

As we indicated above, temperature is the key factor in the non-Markovian system 
coefficients. In Fig. 1, we plot the coherence function A{t) vs "temperature fc^T" vs Uot 
in r = 0.1, and the initial state (xq, yo, Zq) = for the quantum system 

without control effect, 

^p(t) = -t[Hsit),p]+T,it)V[a-]p + T2it)V[a+]p. (26) 

From Fig. 1 we can compare the non-Markovian coherence dynamics with the 
Markovian one clearly. The up figure is the non-Markovian one from which we can see 
the oscillation of the A(t). Moreover, at the low temperature especially temperature 
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r=0.1, non-Markovian coherence evolution 




r=0.1 , Marl^ovian coherence evolution 




Figure 1. (Color online) Comparing the non-Markovian coherence dynamics with 
the Markovian one by the time evolution of K{t) as a function of ksT for initial state 

(xo, 2/0, ^o) = ^, ^) > ^ - 0.1, = 1, and ^ o.Ol. 

the non-Markovian effect is faint, as the temperature rises, the non-Markovian 
becomes more and more obvious, while the Markovian one decays exponentially. This 
phenomenon embodies the non-Markovian effect, which is evidently different from the 
Markovian property. The reason is that due to the non-Markovian memory effect, 
particularly A(t) < in Eq ( 12T|) . the coherence oscillates. With A(t) — 7(t) > the 
quantum system coherence descended whilst A(t) — 7(t) < the coherence ascended. In 
high temperature the Markovian quantum system decays exponentially and vanish only 
asymptotically, but in the non-Markovian system the coherence A(t) oscillates, which 
is evidently different from the Markovian. In this case the non-Markovian property 
becomes evidently. From ( 12T]) . we have learned that in the strong coupling, high 
temperature and high cutoff frequency regimes, the considerable back-action of the 
non- Markovian reservoir effectively counteracts the dissipation. 

4. Example: quantum decoherence control 

As is introduced in the introduction, quantum computing and quantum communication 
have attracted a lot of attention due to their promising applications such as the 
speedup of classical computations and secure key distributions. Although the physical 
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Figure 2. (Color online) Dynamics of coherence function K{t) of non-Markovian 
optimal decoherence control (magenta solid line) , non-Markovian without control (red 
dotted line), and Markovian without control (blue dashed line) for difference cases: 
(a) r = 0.5, ksT = 1, (b)r = 3, fcsT = 1, (c)r = 0.5, hgT = 10, and (d) 
r = 3, ksT ~ 10 respectively. The other parameters are chosen as M = 0.05, and 
77 = 1. 



implementation of basic quantum information processors has been reported recently, the 
realization of powerful and useable devices is still a challenging and yet unresolved task. 
A major difficulty arises from the coupling of a quantum system to its environment that 
leads to decoherence. Various methods have been proposed to reduce this unexpected 
effect in the past decade, which can be divided into coherent and incoherent control, 
according to how the controls enter the dynamics. However, the above control strategies 
render the quantum systems are always neglect the quantum measurement back-action 
or study simple systems with the Markovian approximation. In this paper, we restrict 
our discussion to the no-Markovian open quantum system, and consider its dynamics 
with feedback control and therefore the dynamics obeys the no-Markovian master 
equation (fT2l) . 

To quantify the decoherence dynamics of the qubit, we apply the following quantity 
A(t) which is determined by the off-diagonal elements of the reduced density matrix. 

|pi2(t)|=A(t)|pi2(0)|. (27) 
In the Bloch representation, we have 

V 2;^ + 
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The decoherence factor maintains unity when the reservoir is absent and vanishes for the 
case of completely decoherence. For definiteness, we consider the following initial pure 
state of the qubit |'?/^(0)) = a|+)+/3|— ). To use the optimal control method fll5lll6fl9|20l) 
we set the optimal control target pt as the free evolution px = —^[uoaz^pr]- It is easy 
to solve the equation with Bloch representation xrit) = xocoswot — yo sin uot,yT(t) = 
xq sin uot + yo coscuot, zxit) = zq. 

To demonstrate the effectiveness of our optimal control strategy, we present 
numerical simulations with the initial state (xq, yo, zq) = and coupling 

constant = 0.01, optimal control weighting factor 6 = 1, and = 1 as the norm 
unit. Moreover, we regard the temperature as a key factor in decoherence process. 
Another reservoir parameter playing a key role in the dynamics of the system is the 
ratio r = Uc/ujq between the reservoir cutoff frequency Uc and the system oscillator 
frequency Uq. In Fig. 2, we plot the dynamics of coherence function A(t) of non- 
Markovian optimal decoherence control (magenta solid line), non-Markovian without 
control (red dotted line), and Markovian without control (blue dashed line) for difference 
cases: (a) r = 0.5, ksT = 1, (b)r = 3, ksT = 1, (c)r = 0.5, ksT = 10, and (d) 
r = 3, kBT = 10 respectively. The other parameters are chosen as M = 0.05, and 
7] = 1. From Fig. 2, it is worth noting that as increasing the ratio r, or increasing the 
temperature ksT the coherence lasting time becomes shorter and shorter. From Fig.l 
we have known that non-Markovian reservoir has dual effects on the qubit: dissipation 
and back-action. The dissipation effect exhausts the coherence of the qubit, whereas 
the back-action one revives it. Here we still can observe the dynamical mechanism of 
the non-Markovian effect: the non-Markovian system's coherence lasting time is always 
longer than the Markovian one. From the simulation results, the coherence function 
A{t) will be completely lost in the absence of control neither non-Markovian system 
(dotted line) nor Markovian system (dashed line). However, the feedback control steers 
it to a stationary stochastic process which fluctuates around the target. 

5. Conclusions 

In conclusion, we have investigated the problem of optimal control of non-Markovian 
open quantum system via feedback. At first we analyzed the non-Markovian quantum 
system and the master equation. In general, the reduction of the degrees of freedom 
in the effective description of the open system results in non-Markovian behavior. The 
non-Markovian master equation for the quantum system thus supports to investigation 
of non-Markovian effects beyond the Born-Markovian approximation. In this paper we 
make a thorough examination of the difference between the Markovian system and 
the non-Markovian one. The main difference is that in the non-Markovian master 
equation, one or several of the dissipation coefficients become temporarily negative which 
expresses the presence of strong memory effects in the reduced system dynamics. From 
our analytic and numerical results, we find that the non-Markovian reservoir has dual 
effects on the qubit: dissipation and backaction. The dissipation effect exhausts the 
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coherence of the qubit, whereas the backaction one revives it. In the strong couphng, 
high temperature and high cutoff frequency regimes, the considerable backaction of the 
non-Markovian reservoir effectively counteracts the dissipation. 

Based on the non-Markovian master equation we analyzed the optimal control 
problem via feedback. For the quantum system is typically different to the classical one, 
the quantum measurement changes the dynamical evolution. Hence, we consider the 
quantum weak measurement and the corresponding master equation is the stochastic 
one. Moreover, we designed the control Hamiltonian with the control laws attained by 
the stochastic optimal control problem and the corresponding optimal principle. Usually 
this kind of problem is difficult to be analytically solved. We considered this problem 
in the non-Markovian two-level system. Through transforming its master equation into 
the Bloch vector representation we obtained the corresponding differential equation with 
two-sided boundary values. 

At last, we considered the exact decoherence dynamics of a qubit in a dissipative 
reservoir composed of harmonic oscillators, and demonstrated the effectiveness of our 
optimal control strategy. Obviously, the coherence function will be completely lost in the 
absence of control neither non-Markovian system nor Markovian system. However, the 
feedback control steers it to a stationary stochastic process which fluctuates around the 
target. In this case the decoherence can be controlled effectively, which may indicates 
that the decoherence rate can be slowed down and decoherence time can be delayed 
through design engineered reservoirs. 
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